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X-ray binaries are composed of a normal star in orbit around a neutron star or stellar-mass 
black hole. Radio and X-ray observations have led to the presumption that some X-ray bi- 
naries called microquasars behave as scaled down active galactic nuclei. Microquasars have 
resolved radio emission that is thought to arise from a relativistic outflow akin to active galac- 
tic nuclei jets, in which particles can be accelerated to large energies. Very high energy 7-rays 
produced by the interactions of these particles have been observed from several active galactic 
nuclei. Using the High Energy Stereoscopic System, we find evidence for gamma-ray emission 
>100 GeV from a candidate microquasar, LS 5039, showing that particles are also accelerated 
to very high energies in these systems. 
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High resolution radio maps of X-ray binaries (XRB) have revealed powerful outflows that are similar to the 
relativistic jets seen in active galactic nuclei (AGN) (7]|2j. In both cases, the radio emission is due to synchrotron 
radiation from particles accelerated to high energies. These outflows probably result from the accretion of material 
onto the compact object, albeit on vastly different scales: the mass (size) of black holes in AGN is at least 10 6 
times larger than that of compact objects in XRB. Hence, XRB with resolved radio emission have been dubbed 
microquasars, reflecting the suspicion of some fundamental scaling with compact object mass. 

The kinship should be most evident close to the black hole, where the jet is launched and where the available 
energy reservoir to accelerate particles is largest. In AGN the particles can reach energies such that their non- 
thermal emission extends to the GeV-TeV 7-ray regime, via Compton upscattering of ambient photons or as a 
result of high energy hadron interactions. Because of relativistic bulk motion, this emission is most easily seen in 
blazars, where the AGN jet is aligned close to the line-of-sight. Very high energy (VHE) 7-rays are to be expected 
from some XRB if the physical processes in the vicinity of the compact object are indeed analogous. However, 
previous observations of VHE emission from XRB were inconclusive (0- 

Two XRB have resolved radio emission in the 0.001 arcsecond range, which is presumed associated with a 
relativistic jet, and possible counterparts in the MeV-GeV domain I0|5||6j. LS 5039 (RX J1826. 2-1450) and LSI 
+61°303 (V615 Cas) are each composed of a massive star in an eccentric orbit around an undetermined compact 
object 1012). Their proposed 7-ray counterparts are respectively localised to 0.5° (3EG J1824-1514) and 0.2° 
(3EG J0241+6103) in the Third Energetic Gamma-Ray Experiment Telescope (EGRET) catalogue (9i, so that 
the association cannot be considered firm on the basis of positional coincidence alone. The systems are quite 
inconspicuous in x-rays, with low variability and luminosity of ~ 10 34 ergs -1 at 1 to 10 keV M(M11\ . which is 
about one-tenth their luminosity above 100 MeV, assuming that the EGRET sources are indeed counterparts. The 
7-ray spectra measured by EGRET are hard, with photon indexes T close to 2, which suggests that emission could 
extend to the >100 GeV regime where atmospheric Cherenkov Telescope (ACT) arrays operate. Constraining 
the emission cutoff energy provides important clues as to the physics of the 7-ray source. Furthermore, the XRB 
associations can be rigorously tested by the superior angular resolution of ACTs. 

Located in the Southern Hemisphere, LS 5039 is ideally accessible to the High Energy Stereoscopic System 
(H.E.S.S.). H.E.S.S. is an ACT array of four telescopes located in Namibia, each equipped with a 107 m 2 mirror 
and a 960-photomultiplier tube camera U2\\13\U4l . The telescopes image the Cherenkov light from showers of 
particles created when VHE 7-rays and cosmic rays enter the atmosphere. A central trigger selects only showers 
seen by at least two telescopes. The combination of high resolution imaging and stereoscopic shower reconstruc- 
tion allows efficient rejection of the background from cosmic ray-initiated showers. The H.E.S.S. 5 a sensitivity 
above 100 GeV reaches 1 % of the Crab Nebula flux after 25 hours of observations close to the zenith. The direction 
of each 7-ray shower is determined to 0.1° accuracy, which enables localization of sources in the arc min range 
within the 5° field of view. 

The Galactic Plane survey carried out by the full H.E.S.S. array in the summer of 2004 testifies to the perfor- 
mance of the instrument. In \15\ we reported on the discovery of eight extended VHE 7-ray sources within ±30° 
of the Galactic Center and ±3° of the Plane. After standard quality selection, a total of 25 pointings (10.5 h live 
time) taken during the scan were found to cover the position of LS 5039. The data were independently analysed 
with two separate calibration pipelines (76i and several different reconstruction methods, all of which were in 
excellent agreement with each other. The results presented here are based on a maximum likelihood adjustment of 
a shower model to the observed images to obtain the direction, impact parameter, and energy of the primary il7t . 
The likelihood adjustment also provides for each event a probability that is used to select the 7-ray-like events. 
An image size cut of 60 photoelectrons was applied to avoid systematic effects due to inhomogeneities of the night 
sky background over the field of view, corresponding to an average post-cut spectroscopic energy threshold of 
220 GeV. 

The reconstructed 7-ray map shows an excess 1° southwest of the previously reported hotspot HESS J1825- 
137 (Fig. 1). The excess on each position is calculated by comparing the number of source events, integrated 
over the instrument point spread function (PSF), against the estimated number of background events in the same 
region. We find a significance of more than la for this new source, denoted HESS J1826-148. The source is 
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point-like with a size upper limit of 50" (la) given by a likelihood fit to a Gaussian source profile folded through 
the detector response. This is actually the only point like source discovered in the Galactic scan (17 5> . The best 
position is a(J2000)=18 ,l 26 m 15 s and 6=-14°49'30" (with statistical and systematic uncertainties of ±32" and 30", 
respectively, comparable to the uncertainties of the other survey sources) (75i. The positional accuracy is limited 
by the presence of an extended nearby source and systematics in observations taken at large offsets. On-axis 
observations should improve the positional error to better than 15". 

The 7-ray spectrum was derived from the comparison of reconstructed event energies (in a circle of 6 arc min 
around the source) to the prediction for a given spectral shape {18). The prediction uses energy resolutions and 
system acceptances derived from simulations, taking into account the zenith angle pointing of the array and the 
off-axis angle of the shower in the field of view for each observation. We find an acceptable fit (chance probability 
of 7% of getting a worse fit) to a power-law with a photon index V = 2.12 ± 0.15 (Fig. 2). The low statistics 
currently limit further investigation of more complex spectral shapes. The average integral flux above 250 GeV is 
5.1 x 10~ 12 photons cm -2 s -1 (with statistical and systematic uncertainties of ±0.8 and ±1.3 photons cm~ 2 s _1 , 
respectively), corresponding to a luminosity of ~ 10 33 erg s _1 at 3 kpc lUJ. Errors on the spectral parameters 
correspond to ler confidence interval. The highest energy measurement is at ~ 4 TeV. 

The positions of the supernova remnant G16. 8-1.1 and the pulsar PSR B1822-14, which are both in the error 
box of the EGRET source and are both plausible 7-ray sources, are inconsistent with the position of HESS J 1826- 
148 (Fig. 1). Production of 7-rays from the interaction of cosmic-rays with the interstellar medium is precluded by 
the low H column density at the location of HESS J 1826- 148 relative to its surroundings {19). The radio position 
of LS 5039 is 84" away from the H.E.S.S. position and well within the 3a confidence region (Fig. 1). We verified 
that there are no other radio or X-ray sources compatible with HESS J1826-148 in the NRAO VLA 1.4 GHz Sky 
Survey {20) and in the XMM/Chandra fields analysed by i21t . The observations are not simultaneous with those 
of H.E.S.S. so we cannot formally exclude a long episode of flaring from a blazar. It would be surprising if this 
blazar were not detected in radio, as blazar emission at this wavelength is persistent at levels of ~ 100 mJy, well 
above the survey sensitivity (3 mJy for a point source). The present evidence largely favors the association of 
LS 5039 with HESS J1826-148 and, by extrapolation, with the unidentified EGRET source 3 EG J1824-1514. The 
present H.E.S.S. data are consistent with a constant flux (Fig. SI). Confirmed 7-ray variability correlated with 
other wavebands or a telltale modulation would fully establish the association. 

Several processes can lead to 7-ray emission in LS 5039. The bulk of the luminosity in the system is emitted 
by the 06.5V stellar companion 10 39 erg s _1 ) at an energy kT+ w 3.5 eV (Fig. 2). The binary separation 

varies from 2i?* to 6i?* (where i? Jr w7T0 11 cm), and the radiation density reaches n*«10 14 photons cm~ 3 close to 
the compact object. These stellar photons can be boosted to 7-ray energies by inverse Compton scattering on VHE 
electrons \22\ . With such radiation densities, the energy loss timescale for electrons in the deep Klein-Nishina 
regime, giving a strict upper limit on the radiative timescale, is ~300 s. A short radiative timescale relative to the 
escape timescale from the system (~ 100 s) implies that inverse Compton emission can be very efficient. 

Accelerating electrons to the required energies may be hindered by such rapid losses. Very high energies may 
be easier to reach for protons, which suffer fewer radiation losses. VHE 7-rays may then be emitted via proton- 
proton interactions with the stellar wind. Assuming a stellar wind with a mass loss rate of 10~ 6 M yr _1 and 
a velocity of 1000 km s _1 , the density at 2i?* is ~ 10 10 protons cm" 3 . For such a density the pp interaction 
timescale is ~ 10 5 s. At least 10 2 /10 5 =0.1% of the protons radiate for free streaming particles, implying a total 
kinetic energy of less than 10 38 erg s _1 . Protons may also interact with stellar photons but the threshold is very 
high ~ 10 17 eV. The ^7 timescale is ~ 10 3 s. 

Photons emitted in the H.E.S.S. energy range can interact with the 3.5 eV stellar radiation before leaving the 
system, producing e + e~ pairs. The cross-section maximum a 11 « 1.7 • 10 -25 cm 2 occurs for 7-rays of energy 
« 100 GeV. The opacity « 20 for a photon travelling a distance r «s 10 12 cm (comparable 

to the binary separation). VHE photons emitted close to the compact object are therefore always well inside the 
'7-photosphere' at which r 77 « 1. This initiates an e + e~ pair cascade that redistributes the absorbed radiation to 
lower frequencies. Gamma-rays at energies below 100 GeV suffer little absorption because of the Wien cutoff of 
the stellar spectrum. At higher energies the opacity decreases as 1AE 7 i23t . The VHE spectrum may therefore be 
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hardened relative to its intrinsic shape. 

The absorption of TeV photons in the system can be mitigated. First, the cross-section threshold and ampli- 
tude are angle-dependent, so that 7-rays emitted in a cone pointing away from the companion are not absorbed. 
Scattering of stellar photons in the wind will tend to isotropize the radiation and diminish this effect. Variations are 
expected because the geometry changes with orbital phase. Second, 7-ray emission need not take place close to 
the compact object. Observations of X-ray emission from XRB jets provide evidence for acceleration of electrons 
to TeV energies on parsec scales \24\ . In LS 5039, acceleration at a shock > 1 AU away from the stellar companion 
would happen beyond the 7-photosphere. 

The association of LS 5039 with HESS J 1826- 148 confirms that like some AGN, XRB are able to accelerate 
particles to at least TeV energies. Shocks from colliding ejecta or from jet - interstellar medium interactions are 
natural candidates. Yet the association with an outflow may be questioned in the absence of a direct detection 
of relativistic motion in radio. The relativistic wind of a young pulsar is a conceivable alternative for particle 
injection J25t . The situation would then resemble that in PSR B 1259-63, a system composed of a radio pulsar 
in a much wider 3.4 yr eccentric orbit around a Be star. TeV emission from PSR B 1259-63 was detected with 
H.E.S.S. close to periastron J26t . The higher wind density in LS 5039 probably smears out any radio pulses. 
The resolved radio emission from LS 5039 would be due to particles (cascade pairs) streaming out of the system. 
Further insights into this system may be gained from combined radio and 7-ray observations: Very Long Baseline 
Array maps attain a spatial resolution of a few AU (<5i tantalizingly close to the 7-photosphere. 

References and Notes 

1. I. F. Mirabel, L. F. Rodriguez, Nature 371, 46 (1994). 

2. A. P. Marscher, et al., Nature 417, 625 (2002). 

3. T. C. Weekes, Space Science Reviews 59, 315 (1992). 

4. P. C. Gregory, A. R. Taylor, Nature 272, 704 (1978). 

5. A. R. Taylor, H. T. Kenny, R. E. Spencer, A. Tzioumis, Astrophys. J. 395, 268 (1992). 

6. J. M. Paredes, J. Martf, M. Ribo, M. Massi, Science 288, 2340 (2000). 

7. J. B. Hutchings, D. Crampton, Publ. Astronom. Soc. Pacific 93, 486 (1981). 

8. M. V. McSwain, et al, Astrophys. J. 600, 927 (2004). 

9. R. C. Hartman, et al, Astrophys. J. Suppl. Ser. 123, 79 (1999). 

10. G. F. Bignami, P. A. Caraveo, R. C. Lamb, T. H. Markert, J. A. Paul, Astrophys. J. 247, L85 (1981). 

1 1. C. Motch, F. Haberl, K. Dennerl, M. Pakull, E. Janot-Pacheco, Astron. Astrophys. 323, 853 (1997). 

12. K. Bernlohr, et al, Astroparticle Physics 20, 1 1 1 (2003). 

13. J. A. Hinton, New Astronomy Review 48, 331 (2004). 

14. S. Funk, et al, Astroparticle Physics 22, 285 (2004). 

15. F. Aharonian, et al, Science 307, 1938 (2005). 

16. F. Aharonian, et al, Astroparticle Physics 22, 109 (2004). 

17. M. de Naurois, et al, Proc. of the 28 th ICRC (Tsukuba) (2003), p. 2907. 

18. F. Piron, et al, A&A 374, 895 (2001). 

19. M. Ribo, et al, Astron. Astrophys. 384, 954 (2002). 



5 



20. J. J. Condon, et ai, Astronomical Journal 115, 1693 (1998). 

21. A. Martocchia, C. Motch, I. Negueruela, Astron. Astrophys. 430, 245 (2005). 

22. V. Bosch-Ramon, J. M. Paredes, Astron. Astrophys. 417, 1075 (2004). 

23. R. J. Gould, G. P. Schreder, Phys. Rev. 155, 1408 (1967). 

24. S. Corbel, et ai, Science 298, 196 (2002). 

25. L. Maraschi, A. Treves, Mon. Not. R. Astron. Soc. 194, IP (1981). 

26. M. Beilicke, M. Ouchrif, G. Rowell, S. Schlenker, IAUC 8300, 2 (2004). 

27. N. J. Tasker, J. J. Condon, A. E. Wright, M. R. Griffith, Astronom. J. 107, 21 15 (1994). 

28. M. Ribo, P. Reig, J. Marti, J. M. Paredes, Astron. Astrophys. 347, 518 (1999). 

29. J. Marti, J. M. Paredes, M. Ribo, Astron. Astrophys. 338, L71 (1998). 

30. J. S. Clark, et ai, Astron. Astrophys. 376, 476 (2001). 

3 1 . The support of the Namibian authorities and of the University of Namibia in facilitating the construction and 
operation of H.E.S.S. is gratefully acknowledged, as is the support by the German Ministry for Education 
and Research (BMBF), the Max Planck Society, the French Ministry for Research, the CNRS-IN2P3 and the 
Astroparticle Interdisciplinary Programme of the CNRS, the U.K. Particle Physics and Astronomy Research 
Council (PPARC), the IPNP of Charles University, the South African Department of Science and Technology 
and National Research Foundation, and by the University of Namibia. We appreciate the excellent work of the 
technical support staff in Berlin, Durham, Hamburg, Heidelberg, Palaiseau, Paris, Saclay, and in Namibia in 
the construction and operation of the equipment. 

Supporting Online Material: Fig. S 1 



6 




Figure 1 : Map of excess 7-ray emission in units of counts for the region around LS 5039. The 
map has been smoothed by the point spread function. The white ellipse shows the 3cr confidence 
regionfor HESS Jl 826- 148. The radio emission from theSNRG16. 8-1.1 is represented by gray 
contours (0.05, 0.1, 0.2, 0.3, 0.4 and 0.5 Jy/beam) obtained from the Parkes-MIT-NRAO 6 cm 
radio survey map (12 71 12#l) . The yellow contours show the 68%, 95% and 99% confidence level 
region of the EGRET source 3EG J 1824- 15 14. The green star marks the position of the radio 
source associated with the microquasar LS 5039. HESS Jl 825-137 is discussed in (15). 



7 




Figure 2: Spectral Energy Distribution of LS 5039 including the spectrum of HESS Jl 826- 
148 (points in black, power law fit in red). The average radio, IR, optical and X-ray fluxes 
are shown in blue d2Pll2ffll.ffill271) . Optical fluxes are not dereddened. The two X-ray spectra 
correspond to the historical 1998 high (Rossi X-ray Timing Explorer, RXTE) and 2003 low (X- 
ray Multi-mirror Mission, XMM) flux observations of LS 5039. The multi-year average flux 
above 100 MeV from the EGRET source 3EG J1824-1514 is shown in gray ©. 



8 




53140 53160 53180 53200 53220 

Modified Julian Day 



53240 



53260 



Figure 3: ONLINE MATERIAL. Night by night integral flux above 250 GeV of HESS J1826- 
148. Taking into account a 20% systematic uncertainty, the hypothesis of a constant flux over 
the four months time span is acceptable at the 15% confidence level. The H.E.S.S. data covers 
~ 27 cycles of the 4.4 day orbital period roughly uniformely in phase. No periodic variations 
are apparent when folding the data using the orbital ephemeris of ©. The highest peak in a 
2-100 day Lomb-Scargle periodogram of the run by run fluxes has a significance <2 a. The 
statistics are insufficient to study daily variations in the spectral parameters. 
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